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^ Abstract 

B 

. We have studied the deterministic dynamics of underdamped single and mul- 

O ■ tiparticle ratchets associated with current reversal, as a function of the ampli- 

tude of the external driving force. Two experimentally inspired methods are 
used. In the first method the same initial condition is used for each new value 
of the amplitude. In the second method the last position and velocity is used 
as the new initial condition when the amplitude is changed. The two methods 
are found to be complementary for control of current reversal, because the 
first one elucidates the existence of different attractors and gives information 
about their basins of attraction, while the second method, although history 
dependent, shows the locking process. We show that control of current rever- 
sals in deterministic inertia ratchets is possible as a consequence of a locking 
process associated with different mean velocity attractors. An unlocking effect 
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is produced when a chaos to order transition hmits the control range. 
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Stochastic models known as thermal ratchets or correlation ratchets |T]] , in which a nonzero 
net drift velocity may be obtained from fluctuations interacting with broken symmetry struc- 
tures 0, have recently received much attention. This interest is due to the possible applica- 
tions of these models for understanding molecular motors [Q, nanoscale friction ||^, surface 
smoothening [Q, coupled Josephson junctions 0, optical ratchets and directed motion of 
laser cooled atoms [0 , and mass separation and trapping schemes at the microscale . The 
fluctuations that produce the net transport are usually associated with noise, but they may 
arise also in absence of noise, with additive forcing, in over damped deterministic systems 
0, overdamped quenched systems and in underdamped systems UTIf 



Inertial ratchets, even in the absence of noise have a very complex dynamics, including 



chaotic motion |11,12|. This deterministically induced chaos mimics, to some extent, the role 



of noise [0 , changing, on the other hand, some of the basic properties of thermal ratchets. 



For example, inertial ratchets can exhibit multiple reversals in the current direction |jTT],|T2[. 
The direction depends on the amount of friction and inertia, which makes it especially 
interesting for technological applications such as microscopic particle separation [Q. 

Jung et al ||TT| studied the case of an underdamped particle periodically driven in an 
asymmetric potential without noise and found multiple current reversals varying with the 
intensity of the external perturbation. They characterized the motion by cumulants of the 
contracted, time-dependent solution of the Liouville equation and distinguished regular from 
chaotic transport. 

Several attempts to find the mechanism causing these current inversions have been made. 
Mateos analyzed the relation between the bifurcation diagram and the current. He 



conjectured that the origin of the current reversal is the bifurcation from a chaotic to a 
periodic regime. Close to this bifurcation he observed trajectories revealing intermittent 
chaos and anomalous deterministic diffusion. Barbi et al. related the transport properties 
to phase locking dynamics. They interpreted the current reversals in terms of different 
stability properties of the periodic rotating orbits and reported cases where current reversals 
appear also in the absence of a bifurcation from a chaotic to a periodic motion. Although the 



origin of current inversions seems to be relatively clear, this has not been enough to propose 
a method to control the multiple current reversals which is important technologically. 

The aim of this paper is to elucidate the relationship between bifurcation diagrams, 
phase locked dynamics and transport phenomena for an underdamped deterministic ratchet 
without noise in order to find a way to control the current reversals. 

Specifically we study the dimensionless equation of motion: 

ex + 7X = cos(x) + /i cos(2x) + F sin(c<ji(:). (1) 

Here, e is the mass of the particle, 7 is the damping coefficient, F and uj are respectively 
the amplitude and frequency of an external oscillatory forcing. The asymmetric potential is 
given by 

U{x) = — sin(x) — — sin(2a;). 

Numerical solutions of Eq. (1) are obtained using a variable step Runge-Kutta-Fehlberg 
method [ITg. We fixed e = 1.1009, 7 = 0.1109, /i = 0.5 and u = 0.67 and studied the 



behavior of the system as a function of F. These parameters were chosen because they are 
in the same region as used by Barbi and Salerno [|T^ . 

Let us start analyzing the case of only one particle with a specific initial condition. The 
values of x{t) and v{t) are sampled with a sampling period = T/20 where T is the period 
of the driving force (i.e. T = 2tt/uj). 

The mean velocity of the particle is defined as: 

'( 1 VT ' 

where the number of periods of the transitory was chosen as nfran=400 and the average 
was done until nmax=500 periods. 

There are at least two different ways to do a real experiment to study the effect of 
the variation of the strength of the external force F. One way is to let the particle evolve 
starting from the same initial condition when F changes. Another possible way is to change 



r in the middle of the trajectory of the particle. Simulations were carried out with both 
methods, that will be called method I and II respectively, corresponding to the two possible 
experimental realizations. 

Let us first present the results of method I. The normalized mean velocity < v > /v^j , 
with v^^ — T^/T, where is the spatial period of the potential is shown in figure la for a 
particle with initial condition (a;o=0, ^0=1) as a function of the amplitude of the external 
oscillatory forcing F. The simulation was done by resetting the initial condition to the 
same values (xo=0, vq—1) when F was changed. The corresponding bifurcation diagram is 
plotted in figure lb. Figures Ic and Id are enlargements for F e [1, 1.05]. There are several 
important remarks concerning these figures: a) there exist several inversions va. < v > inside 
regions of the same locking characteristics, which means regions where v is periodic with 
the same period conmensurate with T, (see for example four of such inversions in F between 
1.01 and 1.015); b) there are regions of different locking characteristics and the same mean 
velocity (see for example the region with F between 1.015 and 1.035 and the region near 
F=1.045); c) the values of F where the inversions m. < v > take place are strongly dependent 
on initial conditions. 

To understand these characteristics, the trajectories are shown in figures 2 to 4. In 
figures 2a and 2b we show the case of two trajectories corresponding to different locking 
zones (F=1.015 and F=1.035). Both trajectories (2a and 2b) show the same net transport 
given by a straight line with slope < v >= but the oscillations x over this straight line 
are different as can be seen in figures 2c and 2d where the phase spaces {x/Tx,v/v^^) are 
shown. 

In figure 3 the case of a mean velocity reversal inside the period- 1 zone is shown; the 
complete trajectories for F=1.01 and F=1.011 are drawn in 3a and 3b respectively. In 3c 
and 3d we show the transitory with greater detail. 

Finally in figure 4 the case F=1.527 corresponding to chaotic motion is analyzed. In 
figure 4a we have plotted the trajectory x{t) showing clearly that we are in the presence of a 
net-transport phenomenon and in 4b we show the phase space with the chaotic oscillations 
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superimposed. 

The above results clearly show that with method I the bifurcation diagram gives the 
behavior of the oscillations superimposed on the mean motion and not the net-transport 
movement. The phase locking analysis is a useful tool for the study of the synchronization 
between auto-oscillatory systems with an external periodic driving force []T^. In that case 
the variables x and v are both periodic (i.e. S^) and the net transport is automatically 
discarded. 

Figs. 5 and 6 show the normalized mean velocity and the bifurcation diagram obtained 
with method II, which means taking last position and velocity of the trajectory of the 
previous F as the initial condition for the new F. The step AF is 0.001 and the particle 
evolves during 499T for each F. Fig. 5 was obtained beginning with F=0. 89665 as the initial 
F. The normalized mean velocity remains unchanged equal to 1 until F approximately 1.05. 
Fig. 6 was obtained beginning with F=0. 89666 as the initial F. The normalized mean 
velocity remains unchanged at -1 until F is approximately 1.08. The initial F was chosen 
so that there is a current reversal between them as can be seen in Fig. 1 c). In Figs. 5a) 
and 6a) there is no current reversal in the range [0.9,1.05], while there are many current 
reversals in the same zone of Fig. 1 a). On the other hand, the bifurcation diagrams of Figs. 
5b) and 6b) are continuous but different from each other, furthermore none of them shows 
the abrupt changes appearing in the same zone of Fig. lb). 

These results may be explained if there are two attractors corresponding to mean veloc- 
ities +V;^ and —fa;. In each of the simulations obtained by method II the particle remains 
locked into one of the attractors. The current reversals of Fig. 1 are due to the fact that the 
border between the basins of attraction changes with F. For F in some region around F=l 
the mean velocity may be positive or negative depending on which basin of attraction the 
initial condition belongs to. This sensitivity to the initial condition may be used to control 
current reversals by means of a small perturbation in F when the particle is near the border 
of a basin of attraction. An experiment may be envisioned where a digital generator is used 
as the external force, in that case the step in F and the precise time when the change is 

6 



applied may be tuned to control the trajectory of each particle. As an example in Fig. 7 we 
show the case of a particle starting at a;o=0 and Vq=1 with F = 0.89665. This particle reaches 
a final value of < f > /f^=l as can be seen in Fig. 1. The bold curve in Fig. 7 shows the 
effect of changing F to 0.89675 at time t = 498. 7T , when the particle is at x = 3.1097 + 2n7r 
with an instantaneous velocity v = 0.911662. Due to this change, movement of the particle 
is reversed. On the contrary, if the same change in F is produced at t = 499T, when the 
particle is at x = 1.0597 + 2mT with an instantaneous velocity v = —1.82108, the inversion 
does not occur as is shown with the thin curve. 

These above results were obtained by studying the trajectory of only one particle, but 
ratchet transport is essentially stochastic. This was already pointed out by Feynman [[18 



when he used the ratchet to introduce the Second Law of Thermodynamics. In deterministic 
inertial ratchets deterministically induced chaos mimics the role of noise |jl5| and this calls for 



the use of a time dependent probability measure, as has been previously done in [|TT],|T0|,|T4 

We use a collection of particles with different initial conditions to study the transport 
phenomena. We work with an ensemble of N particles having identical initial velocities 
Vq but initial positions equally distributed in the range [xmin, a;niax]- The initial probability 
density is given by: 

p{x, V,0) = 6{V - Vq) [H{x - Xmin) - H{x - Xmax)] , 

where H{x) is the step function. 

The normalized mean-velocity of the ensemble is: 

1 ^ 

As for the one particle case we first show the results of simulations with method I, which 
in the packet of particles means returning to the initial condition p{x, v, 0) when F changes. 

The results shown in figure 8 were obtained with A^=200 and [a;min,2;max] = [5.08, 11.35] 
and Vo = 0.01 corresponding to a well spread-out packet of particles with initial positions 
between two maxima of the potential. The corresponding bifurcation diagram is also shown. 



There is only one current reversal at F ^ 1.05, where an order- disorder transition takes place. 
About three quarters of the particles in the ensemble have mean velocity < v >j= —v^ 
and the remaining quarter have < v >i= giving a normalized ensemble mean velocity 
<< V » /v^ ~ —0.5. These results seem to agree with Mateos' conjecture that 
current inversion is associated with order- disorder transition in the bifurcation diagram. 

In order to obtain the behavior corresponding to a narrow initial packet, we work with 
the initial condition [xmin,a;max] = [5.08,5.09] and vq = 0.01. In this ensemble particles are 
initially located near the maximum of the potential and having a small initial velocity. Fig. 
9 shows a current inversion around F = 1.015. The corresponding region at the bifurcation 
diagram has no order-disorder transition, contradicting Mateos' conjecture. However, the 
same current reversal, associated with an order-disorder transition, which was obtained for 
the case of a wide packet for F ~ 1.05, does takes place. 

Initial sets of particles with positions near the minimum of the potential and sets of 
particles with identical initial positions and different initial velocities equally distributed 
were also studied with qualitatively similar results. 

It is possible to use method II to control the normalized mean velocity of the packet. 
For example to obtain a normalized mean velocity << V » /v^^ = —0.5 which is the 
minimum normalized mean velocity of Fig. 9 a), with the same narrow packet used above 
we use method II starting with Tmin = 1-013 and increasing it up to F = 1.06 in 100 equal 
steps. The time period for each step in F was conmensurate with T. The normalized mean 
velocity as a function of F is shown in Fig. 10 a). The corresponding bifurcation diagram is 
shown in Fig. 10b). As F changes each particle remains locked to its mean velocity < v >i 
corresponding to Tmin and the packet also remains locked to its initial mean velocity. This 
behavior persists until F = 1.0495 where the mean velocity drops to << V » /v^ = — 1. 
The unlocking effect corresponds to a chaos to order transition. If the simulation starts with 
any Tmin producing a positive velocity locking for the packet the unlocking effect produces 
a current reversal as Mateos found. 

In summary, we studied the variation of the mean velocity of one particle and of a packet 
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of particles with the amphtude of the external force F in two ways. The first one consists 
of returning to the initial conditions whenever F is changed. In this way many current 
reversals appear which are not associated with the bifurcations in the bifurcation diagram. 
The second way is to take as initial condition the last position and velocity of the previous 
trajectory for a given F. When this method is used, the mean velocity remains locked in 
the periodic zones of the bifurcation diagram. The differences between the results obtained 
with both methods may be explained by the presence of at least two attractors associated 
with positive and negative mean velocities. The two methods are complementary as possible 
mechanisms for the control of current reversal, because method I reveals the existence of 
different attractors and gives information about their basins of attraction, while method II, 
although history dependent, shows the locking process. 

We conclude that control of current reversals in deterministic inertial ratchets is possible 
as a consequence of a locking process associated with different mean velocity attractors. An 
unlocking effect is produced when a chaos to order transition limits the control range. 
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FIGURES 

FIG. 1. a) Normalized mean velocity of a particle with initial condition = 0,vo = 1 as a 
function of T. Initial condition is reset whenever F changes (Method I), b) Velocity bifurcation 
diagram for a particle with initial condition xq = 0,uo = 1 as a function of F. Initial condition 
is reset whenever F changes (Method I), c) Enlargement of Fig. la) for F between [1, 1.05] d) 
Enlargement of Fig. lb) for F between [1, 1.05] 

FIG. 2. a) Trajectory of a particle with initial condition xq = 0,uo = 1 for F = 1.015 b) 
Trajectory of a particle with initial condition = 0,vo = 1 for F = 1.035 c) Phase Space 
corresponding to oscillations superimposed on the net motion of a particle with initial condition 
xo = 0,^0 = 1 for F = 1.015 d) Phase Space corresponding to oscillations superimposed on the net 
motion of a particle with initial condition xq = 0, = 1 for F = 1.035 

FIG. 3. a) Trajectory of a particle with initial condition xq = 0, uo = 1 for F = 1.010 b) 
Trajectory of a particle with initial condition xq = 0, = 1 for F = 1.011 c) Enlargement of the 
transitory corresponding to the trajectory of a particle with initial condition xq = 0, = 1 for 
F = 1.010 d) Enlargement of the transitory corresponding to the trajectory of a particle with initial 
condition xq = 0,vo = 1 for F = 1.011 

FIG. 4. a) Trajectory of a particle with initial condition xq = 0,vo = 1 for F = 1.527 b) 
Phase Space corresponding to oscillations superimposed to the net motion of a particle with initial 
condition xq = 0,^0 = 1 for F = 1.527 

FIG. 5. a) Normalized mean velocity of a particle with initial condition xq = 0,t;o = 1 as a 
function of F starting at F = 0.89665. The initial condition for a given F is the last position and 
velocity of the trajectory of the previous F (Method II). b) Velocity bifurcation diagram for a 
particle with initial condition xq = 0, = 1 as a function of F starting at F = 0.89665. The initial 
condition for a given F is the last position and velocity of the trajectory of the previous F (Method 
II). 
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FIG. 6. a) Normalized mean velocity of a particle with initial condition xq = 0,vq = 1 as & 
function of T starting at F = 0.89666. The initial condition for a given F is the last position and 
velocity of the trajectory of the previous F (Method II). b) Velocity bifurcation diagram for a 
particle with initial condition xq = 0, uo = 1 as a function of F starting at F = 0.89666. The initial 
condition for a given F is the last position and velocity of the trajectory of the previous F (Method 
II). 

FIG. 7. Example of control of current direction by selecting the time when F is changed from 
F = 0.89665 to F = 0.89666. The bold curve corresponds to a changing time t = 498.7r. The thin 
curve corresponds to a changing time t = 499T. 

FIG. 8. a) Mean velocity of a set of particles with initial positions between two maxima of 
the potential [5.08, 11.35] and initial velocity vq = 0.01 as a function of F. Initial conditions are 
reset whenever F changes (Method I), b) Velocity bifurcation diagram for a particle of the set of 
particles with initial positions between two maxima of the potential [5.08, 11.35] and initial velocity 
vo = 0.01 as a function of F. Initial conditions are reset whenever F changes (Method I). 

FIG. 9. a) Mean velocity of a narrow set of particles centered at the maximum of the potential 
and initial velocity vq = 0.01 as a function of F. Initial conditions are reset whenever F changes 
(Method I), b) Velocity bifurcation diagram for a particle of the set of particles centered at the 
maximum of the potential and initial velocity vq = 0.01 as a function of F. Initial conditions are 
reset whenever F changes (Method I). 

FIG. 10. a) Example of control of mean velocity for a narrow distribution of particles centered 
at the maximum of the potential and initial velocity vq = 0.01 as a function of F starting at 
F = 1.013. The initial conditions for a given F are the last positions and velocities of the trajectories 
of the previous F (Method II) . b) Velocity bifurcation diagram for a particle of the set of particles 
centered at the maximum of the potential and initial velocity vo = 0.01 as a function of F starting 
at F = 1.013. The initial conditions for a given F are the last positions and velocities of the 
trajectories of the previous F (Method II). 
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